Aug. 5, 1995 

Quantum Group, Bethe Ansatz 
and 

Bloch Electrons in a Magnetic Field 

Yasuhiro Hatsugai^*, Mahito Kohmoto^, and Yong-Shi Wu^ 

^ Department of Applied Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113, 

Japan. 

Institute for Solid State Physics, University of Tokyo, 1-22-1 Roppongi Minato-ku, Tokyo 106, 

Japan. 

^ Department of Physics, University of Utah, Salt Lake City, Utah 84112, U.S.A. 



Abstract 

The wave functions for two dimensional Bloch electrons in a uniform mag- 
netic field at the mid-band points are studied with the help of the algebraic 
structure of the quantum group Uq{sl2). A linear combination of its gener- 
ators gives the Hamiltonian. We obtain analytical and numerical solutions 
for the wave functions by solving the Bethe Ansatz equations, proposed by 
Wiegmann and Zabrodin on the basis of above observation. The semi-classical 
case with the flux per plaquette (j) = 1/Q is analyzed in detail, by exploring 
a structure of the Bethe Ansatz equations. We also reveal the multifractal 
structure of the Bethe Ansatz solutions and corresponding wave functions 

when (p is irrational, such as the golden or silver mean. 
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I. INTRODUCTION 



The effects of a magnetic field on electrons in quantum regime is a fascinating problem 
with very rich physics. The path-dependent geometric phase of the wave function is known 
to play a fundamental role in the problem. For example, the geometrical phase and other 
geometrical considerations, such as the gauge invariance 0] and the effect of edges 0, are 
very crucial to the quantum Hall effect. |^JT[] Moreover, the quantized Hall conductance has 
a topological character, both for a system without boundary (periodic boundary conditions) 
^ and for a system with boundaries (edges) . Thus geometry and topology are of central 
importance in this realistic condensed-matter problem. 

In this paper we will show that a certain exotic algebraic structure plays an important 
role as well, if one includes the effects of a periodic potential in this problem. This is not 
too surprising, since there are two fundamental periods in the problem: One is the original 
period in the geometric phase, and the other is that of the periodic potential. Interplay of 
the two intrinsic periods brings more and interesting structures and, therefore, much richer 
physics into the problem. 

Normally a system with strong periodic potential is modeled by the tight-binding Hamil- 
tonian and the effects of the magnetic field is included by the Pierls substitution. This 
procedure is usually motivated as an approximation. However, the gauge invariance is pre- 
served by the substitution, i.e. both of the original continuum and the lattice systems have 
the same local U{1) symmetry. Thus in our opinions, the Pierls substitution has captured a 
fundamental feature of the problem and, therefore, has a meaning more fundamental than 
as merely an approximation. 

So in the problem of a Bloch electron in a uniform magnetic field, one can use the flux 
per plaquette, 0, to characterized the system. The equations of motion can then be reduced 
to one-dimensional ones, which appear in many different physical contexts, ranging from the 
quantum Hall effect [^,|^ to quasi-periodic systems 0. When is irrational, that is, the 
magnetic flux is incommensurate with the periodic potential, the spectrum is known to have 
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an extremely rich structure like the Cantor set and to exhibit a multifractal behavior. 
Another interesting case is the weak field limit. When the fiux is small, the semi-classical 
treatment of the WKB type is justified. 

Recently a relationship is found between this physical problem and a mathematically new 



concept, the so-called quantum group. Wl\ The essential point is that at least for some spe- 
cial points in the spectrum, e.g. at the mid-band points, the Hamiltonian (in certain gauge) 
can written as a linear combination of generators of a quantum group, called Uq{sl2). On the 
basis of this observation, Wiegmann and Zabrodin has proposed Bethe-Ansatz type equa- 
tions for corresponding wave functions. This opened the possibility of solving the problem 
analytically. Indeed, we have turned the dream into reality: We have found explicit analyt- 



ical solutions to this problem at the first time in the literature. |]T4[ The algebraic structure 
of Uq{sl2) and the associated Bethe Ansatz equations also gives us a new method to solve 
the problem numerically. Especially when we study the irrational limit of a well-organized 
sequence of rational fiuxes, the new method is convenient and beneficial in revealing the 
multifractal behavior, as we briefiy reported also in ref. [Q. 

In this way, it becomes clear that in addition to the underlying geometrical and topo- 
logical structures, the problem of a Bloch electron in magnetic field also has an interesting, 
fundamental algebraic structure hidden in it. Previously there were works in the literature 
about possible relation of the quantum groups to physics. Many of the works, however, deal 
with a sort of artificial (rather realistic) g-deformation of an original physical problem. On 
the other hand, for the Bloch electrons, the quantum group is for a realistic physical prob- 
lem: We are dealing with ordinary electrons, not the ^electrons. Therefore we have shown 
that quantum groups can have a fundamental meaning for observable real phenomena. 

In this paper we are going to report on more details of our published results and to 
present a number of new solutions and results. We will describe our analytical solutions 
in detail and in a self-contained manner. We also present more numerical solutions of the 
Bethe Ansatz equations when the analytical ones are not available. In the next section. 



Sec. 0, we review the Bloch electron and the magnetic translations. In Sec. |TT| and |V 



we summarize necessary stuff about the quantum group Uqi^sh), its representations and the 



functional Bethe Ansatz equations following Wiegmann and Zabrodin. [|TT| Then in Sec. 
|V], we explain in detail our way of analytically solving the Bethe Ansatz equations at a 
particular mid-band point _E = 0. Sec. [V^ is devoted to the analysis of the wavefunction by 
using the explicit solutions obtained in the previous Sec. |V[ We also present a multifractal 
analysis of the explicit wavefunctions for irrational fluxes, such as the golden and silver 



means. Finally in Sec. [VI1| , we present numerical solutions to the Bethe Ansatz equations 



at other mid-band points than E = Q. Some technical details are put into several Appendices 
for the convenience of reader. 

II. BLOCK ELECTRONS AND THE MAGNETIC TRANSLATIONS 

A tight-binding electron in a magnetic field is described by the Hamiltonian 

H = T, + Ty + Tl + Tl (2.1) 
where T^, and Ty are the covariant translation operators defined by 

Tx ^ ] ^m.+l,n^m,n(i rn,n ^ (2.2) 
m,n 



m,n 



where Cm,n is the annihilation operator for an electron at site (m, n). The phase factors 
and 9]^ „ are related to a flux per plaquette 0m,n at (m, n) by 

rot(^,„)e = AJl^^ - A,e^^„ = 27r0^,„, (2.4) 

where the difference operators A^; and operate on a lattice function fm,n as 

^xfm,n fm+l,n fm,n (2'5) 
^yfm,n fm,n+l fm,n- (2'6) 

As mentioned in the introduction, H has a local U{1) gauge symmetry, i.e., it is invariant 
under 



Ci VtiCi (2.7) 
e^^^j Vtic'^^^Vt]^ (2.8) 
= 1, Vj = (m,n). (2.9) 

The covariant translations and Ty do not commute with each other. For a one-particle 
state \^m,n) = cj„n|0) which is localized at the site (m, n), their operation satisfies 

T,T,|^^,„) = e'^^'^--"T,Ty\<iIm,n). (2.10) 

Now let us assume that the magnetic field is uniform, so that we can choose a gauge in 
which both A^6'^„ and Ay9^^ are constants independent of the site {m,n). In this case, 
there is an important symmetry of the Hamiltonian, described by the magnetic translation 
operators and Ty which commute with the Hamiltonian 

[H,t] = [H,fy] = 0. (2.11) 

They are explicitly given by 

Tx = '^m+l.ra'-'m.n^ (2.12) 
m,n 

Ty = ^ '^m,n+l'^m,n^ Xm,n ^ (2.13) 
m,n 

where the phases satisfy the relations 

A,x?;.,„ = A.^^^„ (=A,^;;^„ + 2vr0^,„), (2.15) 
AxX\^ = AyQl^^^ (= A.^^^„-2vr0^,„), (2.16) 
AyX\^ = t^yQl,^^. (2.17) 

In the parenthesis, we have used Eq. (p.4|). We can easily solve them as 

X™,„ = + 2v™0^,„ (2.18) 
X:^,„ = - 27rm<^„,n, (2.19) 



which are gauge dependent. 

In the following, let us assume the magnetic field is uniform and rational, that is. 



P/Q with mutually prime integers P and Q. We can take a "diagonal" gauge [T^JTll as 



C,n = vr0(m + n), = -vr0(m + n + 1), (2.20) 

which is derived from the Landau gauge in (1, 1) direction. Correspondingly, we define 
magnetic translations along (1, 1) and (—1,1) directions T+ and T_ respectively as 

TV = f^fy, (2.21) 
f_ = f-^fy. (2.22) 

They commute with the Hamiltonian, [if, T+] = T_] = 0. They do not commute with 
each other. However, one can explicitly check that and T_ commute with each other. 
Thus we can take simultaneous eigenstates of H and and T_, which are specified by the 
momentum in the magnetic Brillouin zone (the Bloch theorem) as 

f^|$(p+,p_)) = e^2QP+|$(p^,p_)), (2.23) 
r_|$(p+,p_)) = e^Pi<l>(p+,p_)), (2.24) 

where |$) is a one particle state spanned as 

MV+,V-)) = E ^™,n(p+,P-)4,„|0). (2.25) 

m,n 

From Eqs. ( |2.23|) and (p.24 we have. 



*m,n(p) = e^^'-("^-")+^P+("^+'"V™+n(p) (2-26) 

tPk+2Qip) = Mp), k = l,---,2Q. (2.27) 
The Schrodinger equation, if|$(p)) = E{p)\^{p)) is written by ipkip) as 

^ ^ ^^^^^^ ^2.28) 



where q is defined as 



q = e'^^ = e^'Q. 



(2.29) 



Using base = [V'i(p); ' ■ ■ ? ^2q(p)]*, Eq. ( |2.28| ) is written in the matrix form as 

H{p)'4>{p) = E{p)%l}{p), with 



(2.30) 



where 2Q x 2Q matrices X and y are given by 



1 
1 



X 



They satisfy relations 
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,F = diag(g^g^■■■,g2«) 



(2.31) 



2Q5 



gXr = YX, X^ = X-\ Y^ = Y 



(2.32) 
(2.33) 



The group generated by X and Y with the relations above is the Heisenberg-Weyl group. 

The solvability of the problem of the Bloch electrons by the Bethe Ansatz lies in the 
fact that the Hamiltonian H{p) has a higher symmetry than the Heisenberg-Weyl group at 
special momentum points. |Tl|] Actually H{p) is written as 



Hmb{P-) = H{p^ = 7r/2,p_) = i{q - q-'){e-'P- pc{B) + e'^-pdC)) 



(2.34) 



at special momentum lines p+ = 7r/2 mod n/Q, usually called as the mid-band lines. (p_ 
is a free parameter, but later we will see that the energy is independent of it.). Here 



PC '■ Uq{sl2) — > M2Q (2(5-dimensional matrices) is a cyclic representation of the quantum 



group Uq{sl2) = {A, B, C, D}. It is explicitly given by 

pc{A) = (2.35) 

Pc{B) = -{q - q-Y'{X - X-')Y, (2.36) 

Pc(C) = -(g - q-Y'y-\X - (2.37) 

pc{D) = g(^-i)/2x (2.38) 



The representation pc will be used to reveal the quantum group structure of the problem 
in the next section. 

III. THE QUANTUM GROUP AND ITS REPRESENTATION 

Historically quantum groups arise from mathematical structures appearing in the exactly 
solvable one-dimensional quantum many-body models. However, the simplest case, Uq{sl2), 
can be understood as an a priori deformation of the ordinary Lie algebra sl2'- Uq{sl2) is ob- 
tained by the g-deformation (i.e. an appropriate "extension" with a new complex parameter 
g, which preserves many good properties) of 5/2 and, conversely, 5^2 is recovered from Uq{sl2) 
by taking a classical limit q^ 1. 

More precisely, Uq{sl2) is an algebra generated by four generators {A^B^C^D} with 
multiplication and addition (with complex coefficients), subject to the following defining 
relations 

(3.1) 

qDB, (3.2) 
qAC, (3.3) 

(3.4) 

This is the g-deformation of the Lie algebra s/2 = {S+iS-^S^}, which is defined by the 
relations [S+, 5_] = 2S'3, [^3, S±] = ±S±. 
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AD^DA^ 1, 

AB = qBA, BD = 

DC^qCD, CA^ 
A^-D^ 



7-1 



It is well known that 5/2 has a {2j + 1) dimensional (spin-j') representation which is 
realized by the differential operators 

PhiSs) ^ zd^ - j, (3.5) 
Ph{S+) ^ z{2j - zd,), (3.6) 
Ph{S-) = d,. (3.7) 

acting on a function "^{z). This finite-dimensional representation has both a highest weight 
and a lowest weight, since ph{S+)z'^^ = ph{S-)z^ = 0; thus, it is a 2j + 1-dimensional 
representation with ^{z) a polynomial of degree 2j. Correspondingly, Uq{sl2) also has a 
highest-weight representation, realized by the g-difference operators: 

Ph{A)^{z)^q-^^{qz), (3.8) 

Ph{B)^{z) = {q- q-Yh[q'^^{q-h) - q-''^{qz)], (3.9) 

PhiCMz) = -{q - q-yz-'[^(q-'z) - ^(qz)], (3.10) 

Ph{D)^if{z)^q^^{q-'z). (3.11) 

It has also a highest weight and a lowest weight: ph{B)z^^ — ph{C)z^ = 0. In the standard 
bases, *(z) = ( n = 0, 1, • • • , Q - 1), they operate as Ph{A)z'' = Ph{D)z'' = 

qj-n^n^ p^i^Q^^n _ _ ^j^^:""*"^, and ph{C)z"' = [n]gz"'~^ where [n] is a g-integer defined as 

N. - ff^, (3-12) 
which approaches to n in the limit 5 — > 1. Thus they can be represented as matrices, given 

by 

Ph{A): diag(g-^g^-^•••,g^-^g^■) (3.13) 
pn{D) : diag(?^ q-'+^ , • • • , q-^+\ q'^) (3.14) 
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Ph{B) : 



Ph{C) : 





[2j-l], 


[1], 

[2], 



[1]. 



(3.15) 







[2j], 




(3.16) 



Notice that the dimension of the matrix above is Q = 2j + 1 and the matrices are all 
tri-diagonal. 

Up to this point, Uq{sl2) is quite analogous to the 5/2- However, when q is a root of 
unity, there is another representation, called as a cyclic representation, whose existence is 
characteristic to Ug{sl2)- (s/2 does not have cyclic representations.) Our case, g^*^ = 1, is 
such a case. Here a 2(5-dimensional cyclic representation is derived from the highest-weight 
representation of spin-j = (Q— 1)/2 as follows. First define a cychc basis {^n: n — 1, - ■ ■ ,2Q 
(mod 2(5)} for the cychc representation (satisfying '^n+2Q = ^n) by 



^ r 



(3.17) 



where '^{z) is the polynomial in the spin-j = {Q — l)/2 highest- weight representation. In 
this realization, operations of the cyclic representation are explicitly given by 



Pc(^)*n-g-(^-')/'*n+l, 

pc{B)^n = -{q - q-Y\q''+'^n+i - q'^'^n-i), 

Pc{C)^n ={q- r')-'(r"*n+l - r"*n-l), 



(3.18) 
(3.19) 
(3.20) 
(3.21) 
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If we write them in matrix form, we get Eqs. ( |2.35| )- ( p.3^ ). We note that pc{B)'^2Q returns 
to a hnear combination of and ^2Q-i- So there is no highest-weight state (or vector) \l/ 
that can be annihilated by pc{B): pc{B)'^ = 0. Similar things happen to pc{C): There is 
no lowest-weight state annihilated by pc{C). These are the reasons why this representation 
is called cyclic. 

Thus, the highest-weight representation of dimension Q is related to the cyclic represen- 
tation of dimension 2Q. This establishes the following results about the spectrum and the 
wavefunction of the Bloch electron at the mid-band momenta. The spectrum is given by 
the eigenvalues of the following Q x Q tri-diagonal hermite matrix. 

[l],e-'P- 

[2]ge-'P- O 
[2],e+'^P- 

[Q-2]ge-'P- 
O [Q- 2]ge+'P- 

[Q - l],e+^^ 



[1] 



[Q - l]ge-''P 


(3.22) 



where we have used the fact [Q — j]g = [j]q. When we take large Q, the eigenvalues of 
the tri-diagonal matrix gives the so-called Hofstadter's butterfly. In Figs. ^ we show eigen 
values of the -^1"^ as a function of = ;^Q'logg for Q = 23 and Q = 401. Since the eigen 
values of the tri-diagonal matrix does not depend on the phase of the off-diagonal matrix 
elements, the energy does not depend on p_. The wavefunction of the Bloch electron at the 
mid band momentum of the j-th band is obtained from the j-th eigen vector {f^}, {m = 
0, 1, ■ ■ ■ , Q — 1) as 



(3.23) 



m=0 



The dimension of the original Hamiltonian Hmb is 2Q and the tri-diagonal functional Hamil- 
tonian is Hj^/^ is Q. Actually, Hmb has doubly degenerate eigen values. The above state 
eq. (|3.23|) is degenerate with the other state 
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*:^ =^(9-")=E (3-24) 

m=0 

These two states have the same energy and they are orthogonal each other as 

2Q-1 1 _ 2Q(m+m') 

E^^n^f= E ^-<'T^7Hwr = o (3-25) 

n=0 0<m,m'<Q~l ^ 

Thus we can obtain all solutions of the original problem from the functional equation. 



IV. THE FUNCTIONAL BETHE ANSATZ FOR THE BLOCH ELECTRONS 

In this section, we derive the functional Bethe Ansatz equation following Wiegmann and 
Zabrodin [|ri|]. We denote p_ by p, in the following. The Schrodinger equation at a mid 
band is 

HMB{p)i^n = i{q - q-')[e-'''pc{B) + e'^pc{C)]^n = ^^n- (4.1) 

Recalling that the cyclic representation is derived from the highest weight representation, 
let us define the functional equation 

i{q-q-')\e-'''pH{B) + e'^Ph{C)]^{z) = 
i{e-'Pqz + e'Pz-^)^{qz) - ~i{q~'e-'Pz + e'^ z~')^ (q-' z) = E^{z). (4.2) 

Once we have a solution of the functional equation Eq. (|4.2| ), we can recover all the solutions 
of Eq. (U) using Eq. {^1^ and Eq. { ^1^ . 



Since "^{z) is a polynomial of finite degree, it can be factorized as 

<i/iz)=Y[iz-zUp))- (4.3) 

m=l 

First let us absorb the p dependence of Eq. ([4.2|) by the gauge transformation 

z = ze-'P, (4.4) 

Zm = Zm{p)e-'P, (4.5) 
Q-1 

^{z) = n - 

m=l 

= (4.6) 

m=l 
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Thus we have 



i{qz + z'^)^!{qz) i{q~^z + z~^)^!{q~^z) = E-^{z). (4.7) 



From Eq. (I^) and Eq. (gj), we get 

+ --')^pf^ - -te-'^ + ^-') "g'^^:7' = (4.8) 

Comparing both sides and use the pole-free condition, then one obtains the Bethe Ansatz 
equations 

|^ = -nf^. '=i.-.Q-i (4.9) 

Q^l + J- m=l ~ l^rn 

Q-l 

E = -i{q-q-^)Y.Zm. (4.10) 

m=l 

From this equation, one can easily see that the energy E is independent of Note that the 
functional solution corresponds to \E'^ is \1'(2;"^), which is not a polynomial and has poles. 



V. EXPLICIT SOLUTIONS OF THE BETHE ANSATZ EQUATIONS FOR S = 

Let us consider from Eq. ([4.8|) . For the zero energy E = Wiegmann and Zabrodin 
showed that '^{z) is given by the so-called continuous q-ultraspherical polynomial [|I5| as 



and (a; g)^ = O^Vl ~ aq^')- [0 The q-ultraspherical polynomials have several interesting 
properties. However, we treat the function '^{z) directly. 

In order to understand the properties of the wavefunctions at = 0, the center of the 
spectrum, we solve the Bethe Ansatz equation ( ^^) explicitly. For E = Eq. ( [4. 81) is 
written as and rewrite it as 

q{z^ + q'^)'^{qz) = q~\z^ + q)'^{q-^z). (5.1) 

First put z = iq^^^'^, then one obtains q{—q + q^^)'^{iq^^'^) = and '^(iq^^'^) = 0. So ig^/^ is 
one of the roots. Next put z = ig"*"^/^, then one obtains g(— g^-|-g~"^)\l'(ig''/^) oc \I'(ig^/^) = 0. 
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Thus ic^l"^ is also a root. We can repeat this procedure and get a series of the roots as ig2m-i/2 
where m = 1, ■ ■ ■ , (Q — l)/2. The restriction of m arises due to the fact that a pre-factor of 
the vanishes at the last step. 

Similarly, put z = iq'^^"^, one obtains another sequence of the roots as ^g-2m+i/2 -^j^gj^g 
m = 1, ■ ■ ■ , {Q — l)/2. Comparing the number of the independent roots with the degree of 
the polynomial, Q — 1, we have obtained a complete set of roots as 

{• 2m-l/2 ^ ^gi27r<^{m-l/4) 
m= l,---,(g-l)/2. (5.2) 
^^-2771+1/2 ^ ^gi27r0(m-l/4) 

They are all on the unit circle. Thus let us write them as Zm{p) = e*^™*^^-* and Zm = e*^™ and 
consider the distribution of 6m by p{6) = limg^oo Q^O, where is a difference between 
adjacent 6m- Since Zm = Zm^^-'', P dependence of the roots are simple shift of 9m{p) = Om + P- 

The restriction of m {1 < m < {Q — l)/2) causes a nontrivial distribution of the roots. 
The phase factor 6m corresponds to a quasi-momentum in the Bethe-Ansatz-solvable models 
in which Q corresponds to the system size. In the following, we consider Q —>■ oo limit which 
corresponds to the thermodynamic limit. 

Consider first the case P = 1. The roots {zm} distribute uniformly on the unit circle 
except near z = ±i. The roots for Q = 89 are shown in Fig. ^. In the semi-classical limit 
Q ^ oo, that is, q ^ I, the distribution function p{9) = limg^oo QA9 is smooth (constant). 
A continuous behavior of p{6) is usually obtained in the Bethe-Ansatz-solvable models in 
which p{6) is determined by an integral equation. 

When the flux is irrational, the situation is quite different. First let us take the flux 
((, = 1/t = - l)/2, where T is the golden mean (-\/5 + l)/2. To reach this flux, we 
consider a sequence of rational fluxes (pk = Pk/Qk, where Qk = P-ik+i, Pk = F^k and is 
a Fibonacci number deflned by F^+i = Fk + Fk-i, Fi = 1, Fq = 1. In this sequenses, Pk 
and Qk are all odd. The two types of roots in Eq. (|5.2| ) are nested and has a complex 
distribution. To gain an insight into the distribution of roots, it is helpful to consider 
the pseudo roots which are deflned also by Eq. ( ^.2|) but with the range of m modifled to 
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m = — (Q — l)/2, ■ ■ ■ , —1, 0. In Figs. ^ we show the distributions of the roots (black points) 
and pseudo roots (gray points) for several 0^, ( k = 1,2,3,4,5,6). Here the radius of the 
unit circle has been scaled so as to show all the cases at once. These figures clearly show 
that there is a branching rule for the true roots (denoted by A) and pseudo roots (by B) as 
follows 

A^ A^B^A^ 

B^ B^A^B^A^B^ 

B^ B'^A^B^ (5.3) 

The initial condition is B^A^B^A^ (cyclic). In the k-th stage of the sequence, the number of 
clusters of the true roots A'^ and are Qk-i — 1 and Pk-i + 1 respectively. This branching 
rule gives rise to a self-similar structure for the distribution p{6) in the limit k oo. 
To characterize the distribution, let us define the generation of roots. According to the 
branching rule ( ^.3] ), each true (pseudo) root branches into a cluster of 3 new true (pseudo) 
roots, each of which in a sense has a parent. At the same time, between these clusters of new 
pseudo (true) roots, there is a pair of new-born true (pseudo) roots which have no parent. 
We assign the generation number to a root so that it is 1 when the root does not have 
a parent, otherwise it is one plus the generation number of its parent. Let us denote the 
number of true (pseudo) roots in the k-th stage with generation g by nA{g, k) {nB{g,k)) {g = 

1, • • ■ , /c). (Then in the special case P = 1, UA^g, k) = UB^g, k) = 6kg-) In the present case, 
there is a recursion formula by the branching rule as nA{g,k) = 3nA{g — l,k — 1), g = 

2, ---,fc, nA{l,k) = 2{Pk-i + 1). Thus nA{g,k) = 2 ■ S^-^P^.g + 1) and nB{g,k) = 
2-35-i(Pfe_,-l). 

Next let us consider the flux (j) = 1/a = \f2 — 1, where a is the so-called silver mean 
+ 1 . To realize this flux in the large Q limit, we consider a sequence of rational fluxes 
0fc = Pk/Qk, where Qk = Gfc+i, Pk = Gk and Gk is defined by Gk+i = 2Gk + Gk-i, G2 = 
1, Gi = 1. In this case, we can apply our Bethe Ansatz technique to each step of the 
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sequence since all Qk and Pk are odd. In Figs. ^ we show the distributions of the roots 
(black points) and pseudo roots (gray points) for several 0^, ( /c = 1, 2, 3, 4, 5, 6). The radius 
of the unit circle are also scaled as before. One can observe a clear branching rule branching. 
It is a little different from the golden mean case and is given by 

BA^BA^B 
BA^BA^BA^B 
AB^'AB^A 

AB^AB^AB^A (5.4) 

The initial condition is AB^AB^ (cyclic). From this branching rule, we have a self similar 
behavior at every two stages. At the odd stage, all the roots and the pseudo roots appear 
as A, A"^, B and B^ and at the even stage, they appear as A^, A^, B^, and B'^. For both 
cases, the distribution function of the roots has self similarity and reflects the self similar 
structure of the original problem. 

The above considerations clearly exemplify the difference of the distributions between the 
semi-classical limit and the incommensurate cases. The distribution of the roots has a self- 
similar structure and the function p{6) is nowhere differentiable in the incommensurate limit, 
while that for the semi-classical limit is smooth. This is characteristic to the incommensurate 
case. 

Another way to characterize the distribution is to map to the dual (reciprocal) space. 
This can be done for arbitrary P and Q analytically. We lift the 6^ to the real axis period- 
ically. On the real axis, the true and pseudo roots occupy lattice points {j/2Q \ j: integer} 
with spacing 1/2Q. Thus we define the Fourier transform of the so-called defining function 
by 

oo 

SQik) = Y: e^'^SQij). (5.5) 

j=-oo 

where Sgi^j) = 1 if there is a true root at j/2Q, otherwise Sq^J) = 0. As shown in Appendix 
A, 
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A^A^ 
A^ ^A^ 

B^B^ 
B'^^B^ 



SQik) = nE^Q ^(^ - fcr), (0 < < 27r), 

V r-=0 



(5.6) 



where 



and 



(g-1) r = 



cos{Pkr/2) 



r = l,---Q-l, 



(5.7) 



(5.8) 



which converges to a usual continuous momentum A; e [0, 27r] in the Q — > oo hmit. In this 
hmit, we have 

1 



\SQ{k)\' 



-, (0 < A; < 2n) 



(5.9) 



4cos(PA;/2)^ 

In the semi-classical limit P — 1 and Q ^ oo, |5'q(A;)P is well defined and behaves smoothly. 
In the incommensurate limit, |S'q(/c)P is a singular function and is not even difi^erentiable 
hence P is infinite. Also it can be shown that the original defining function is given by 



,os(7rPn/g)- 



(5.10) 



VI. EXPLICIT WAVE FUNCTIONS BY THE BETHE ANSATZ SOLUTIONS 

Let us consider wavefunctions which are obained from the explicit solutions of the Bethe 
Ansatz equations. The wave function at site j is given by — "^{q^) and it is written in a 
compact factorized form as 

(Q-l)/2 

= H (5^_/2m-l/))(^,_^-i(2m-l/)) 

m=l 

= {-q^i^q-'^'^'■, q')iQ-mi^q-'''^'; r')(Q-i)/2, (6.1) 



(Q-l)/2 (Q-l)/2 

E E 


(g - 1)/2 




(Q - l)/2 




(6.2) 


i>=0 v'=0 


1/ 


q 


u' 


q 
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where we have used a g-binomial theorem 



11(1 + 9 



N-l-2n, 



n=l 



N 

E 

u=0 



N 



V 



\v\\\li-v\ 



(6.3) 



(6.4) 



It is convenient to shift the site, j = j + Ji, by an amount J\ where J\ is determined by 
PJi = (g - P)/2 (mod 2Q ). Then = at j = 2m, -2m + 1 (m = 1, ■ ■ ■ , (Q - l)/2), 
and the wavefunction is nonzero only at j = 1, 3, ■ ■ ■ , Q and Q + 1, Q + 3, ■ ■ ■ 2Q as shown 
in Appendix B. Thus one has 



1^^ 



Q 



(6.5) 
(6.6) 



[J-I].!!' 

for J = 3, 5, ■ ■ ■ , Q and = l^p where [j],!! = [j] Jj - 2], • ■ ■ [2]„ or[l],. 

When P = 1, we can proceed further analytically. As shown in the Appendix B, we get 
a compact form for the wavefunction in the semi-classical limit, ( P = 1 and Q oo), 

2 



(0 < X < 1), 



(6.7) 



sml^vrxj 

up to a constant factor, where ipi^x) = "^j, x = (2j — 1)/2Q. The squared amplitude of 
the wavefunction is given by the inverse chord distance. The recursion relation = 
|\E'j„2p sin^(7r(j — 2)/Q)/ sin^(7r(j — 1)/Q) obtained from Eq. ( |6.6| ) has played a key role. 
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We can also calculate the finite-size correction near the edge a; ~ or a; ~ 1 for Q ^ oo. 
As shown in the appendix C, it is given by 

2 



\^{x,i+,)\' = C{l)- 

sm 7rx2i+i 

k=l 



(6.8) 
(6.9) 

(6.10) 



where C(0) = 7r/4 = 0.78539..., C(l) = Svr/IG = 0.98174..., C(2) = 8l7r/256 = 
0.99402..., Since sin;z = zn^i(l - n^i(l - 1/4A;2) = n/2. So the finite-size 

correction factor C{1) converges to unity very rapidly and the Eq. ( |6.7| ) is quite accurate 
even for small Q. The norm of the wavefunction is logQ + const, and unnormalizable 
which is characteristic in a critical wavef unctions. In Fig. ^ the amplitudes of the analytic 
wavef unctions, normalized by the peak height, are shown for several values of Q. 

Next let us discuss the case with golden-mean fiux. We plot the analytic results Eq. ( |6.6|) 
in Fig. § for a sequence of rational fiuxes converging to 1/r. One can easily recognize the 
self-similar behavior of the wave function. Each peak branches into three peaks in the 
next stage. Presumably these are the refiection of the self-similar distribution of the roots, 
Eq. (P). 

We have performed a multifractal analysis to investigate the nature of the wavefunc- 
tions. This is useful to distinguish critical wavefunctions from extended wavefunctions. The 
results for critical wavefnctions reveals the multifractal properties. Let us consider a fc-th 
generation Q = Qk = -^sfc+i- First define a probability measure pj {J2jPj = 1 ) of the wave 
function as 

p.-ijn' (6-11) 

^fc = El^/- (6-12) 

i=i 

Next define a Lebesgue measure Ik of each site as 

h = ^. (6.13) 

From Pj and Ik, the singularity of the probability measure is represented by an exponent aj 

as 

P. = C- (6.14) 

The distribution Qk{<^) characterizes the wave function, where f2fc(a)(ia is the number of 
sites whose value of aj lies in the interval [a, a + da]. Since Qk increases exponentially 
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as k increases, so does Vtk- Thus it is natural to introduce the entropy function S{a) to 



characterize the wave function as pO| 



or 



S{a) = hm Sk{a) = hm — logQk{a), (6.15) 

fc— +00 fc— >oo K 



/(«) = (6.16) 



where 

e = hm efc = — hm — log/^- (6-17) 

fc— >oo fc— >oo 

Thus we have ^^(a;) ~ ^k^^"^ ■ calculate /(a) let us define the partition function Zk{r) 
as 

^fcW = Ei>,^ = E/r'"^". (6.18) 

One obtain /(a) = lim^^oo fk{(^) from 

Gfc(r) = ^logZfc(r), (6.19) 

a = --^Gk{r), (6.20) 
Cfc ar 

A(a) = -G'fc(r)+ra. (6.21) 

When the wave function is extended, pj ~ Ik and /(a = 1) at a = 1. On the other hand, 
when it is localized, /(a) consists of two points at a = and a = 00. When /(a) is a 
continuous function, the wave function is critical and has a singular probability measure 
which is a characteristic of the multifractality. We have performed the above multifractal 
analysis. Here we stress the importance of the finite size effects. In order to have a reliable 
/(a), one must perform extraporations of the finite size data, (/(a) obtained from a finite 
system is different from the true /(a) [f7|,pT|.) We have done such calculations. The results 
is shown in Fig. |^ for a golden mean flux. It gives a smooth f{a). This clearly shows 
that this wavefunction is multifractal and critical. We note the striking resemblance of 
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these wavefunctions to that of the Id quasicrystal Fibonacci lattice at the center of the 
spectrum. |18[ The latter was obtained exactly by a different technique and f{a) is obtained 
analytically. [^] 

VII. NUMERICAL SOLUTIONS OF THE BETHE ANSATZ EQUATIONS FOR 

For E ^ cases, we have not being able to obtain analytical results. Thus we investigate 
them. However, the Bethe Ansatz equations are high degree algebraic equations of many 
variables and it is extremely difficult to handle directly even numerically. Thus we use the 
information of the quantum group to reduce the difficulty. Instead solving the Bethe Ansatz 
equations, Eq. (|4.9|) , we construct a polynomial "^^{z) for a mid band energy when \E'-' is 
written as When one write the "^{z) as 

^'iz) = j:vLz"^, (7.1) 

m=0 

[vq, ■ ■ ■ , Vq_i]'^ is given by the j-th eigen vector of the Q x Q tri-diagonal matrix H^^b- The 
roots of the Bethe Ansatz equations are given by the roots of this one variable equation. 
Since the H^mb is real symmetric, '^^{z) is a real polynomial. Thus the roots can be obtained 
by the traditional numerical technique. 

Let us consider the cases P = 1 first. We have calculated the roots of the Bethe Ansatz 
for substantially many different Q. All roots are on the unit circle. It seems that all roots 
of the Bethe Ansatz equations are on the unit circle when cf) = 1/Q with odd Q. In Fig. p, 
we present results for several cases with Q = 61. In Fig. all root for Q = 41 are shown. 
At the highest energy band, the roots are on the right semi-circle of the unit circle and it is 
almost uniformly arranged (not exactly for finite Q). At the second energy band, one root 
appears in the left semi-circle (at z = —1). At the third one, one more roots appear in the 
left. In this way, one of the roots are shifted from right semi-circle to the left semi-circle 
as the band energy decreases. This behavior is clearly observed in Fig. |^. We have also 
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calculated the distribution functions which are plotted in Fig. |T0[ For a finite Q case, we 
define an approximate distribution function pf (6') as 

pQ^e)=Y.5^{o-e^), (7.2) 

m=l 

U^) = --^ L ■ . ), (7.3) 

TT 2 sin(|) — 

where e is a width of the approximate delta functions and we take as e = 2/(5. The distance 
is measured by the chord distance. In Fig. (p!0D , we plot the results for several cases. From 
the results for the finite Q, we can speculate a behavior of the distributions in the Q — >■ oo 
limit. It is likely to have singularities for 7^ states at z = ±z which was proved for = 
state. 
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APPENDIX A: DERIVATION OF THE STRUCTURE FACTOR OF THE BETHE 

ANSATZ ROOTS 

Here we derive Eq. ( p.6D -Eq. ( p.8| ). Using Eq. ( [5. 2] ), one obtains 

00 

j=~oo 
00 0-1 

^ ^ ^ |-gifc{20/+2Pm+i(0-lP)} _^ gjfc{2Qi-2Pm+(0+P)/2}j 



/=— 00 m=l 



00 f 1 _ jPfc(Q-l) 1 _ iPfc(Q-l) >| 

V P*2fc0« J ^ ^ i(Q+3P)fc/2 , ^ ^ i(Q-3P)k/2 I 

I 1 _ gj2Pfc ^ ]^ _ g-i2Pfc ( 



l=—oo 
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The structure factor is a sum of the delta function but the amphtude of the peak has a 
nontrivial feature. At k = we have 



sin Pk 2 



Q-1 r' = 



cos TT-Qf^ 



r' = 2r^0 (A2) 



otherwise 

This leads to Eq. ( ^.GD -Eq. (|5.8| ). In Eq. (|5.7| ), cosPk/2 is smooth function of k E {—11,71] 
when P is finite. It, however, behaves wildly when P ^ 00. This case corresponds to an 
irrational flux. 

APPENDIX B: EXPLICIT WAVE FUNCTIONS IN A CLOSED FORM 

From the explicit solution of the Bethe Ansatz equation, the wave function at site j is 
written as 

{Q-i)/2 

m=l 

{Q-l)/2 
m=l 

= {-q)-K^q"'-': q\Q-l)/2{^q-"'-'■. q-^Q-m- (Bl) 

The two sequences of the roots are obtained explicitly and we have 

(Q-i)/2 1 
^^. = g(Q-i)i J] (i_exp[^27r— {P(2m-j) + (Q-P)/2}]) 

m=l 

X (1 - exp[z27r^{P(-2m - j) + (Q + P)/2}]) (B2) 



Let us define Ji by 



PJi = i(g - P) (mod 2g), (B3) 



then P(Ji + 1) = + P) (mod 2Q) and = 0, for j = 2m + Ji, -2m + Ji + 1. 
(m = 1, ■ ■ ■ (Q — l)/2). Thus we have 
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% = 0, for j = 2m, or - 2m + l(m = 1, • • • (Q - l)/2) (B4) 
where j — j — Ji- 

Now let us consider the amphtude of the wave function. Apparently period of is 

2Q. However, it is not the smallest period. Actually the smallest period is Q because 

{zm} — {^m}- Thus define I hy j — 21 — 1 and consider — ■0j=2i-i — Ji + 21 — 1, 

I — 1,2, - • • ,{Q + l)/2). Then the squared amplitude is 

(Q-i)/2 p p 

m)\'= n |l-exp[.27r{ — (2m -!)--(/ 

X |1 - exp[i27r{^(-2m) - ^{l - 1)}]|^ 

(<3-i)/2 

^ JJ (l_^-2(i-l)^2m-l)(i_^2(/-l)^-(2m-l)) 
m=l 

X (1 - - g2(i-l)^2m). (B5) 

(B6) 

Let us first consider ior I — 1. Using 

,Si = {2m - 1 |m = 1, • • • , (Q - l)/2} = {1, 3, • • • , Q - 4, Q - 2} 

^2 = {-(2m - 1) |m = 1, • • • , (g - l)/2} = {-1, -3, • • • , -(Q - 4), -(Q - 2)} 

= {Q + 2,Q + 4---,2Q-3,2g- 1} (mod 2g) 
^3 = {-2m |m = 1, ■ ■ ■ , (Q - l)/2} = {-2, -4, • • • , -(g - 3), -(g - 1)} 

= {g + 1, g + 3 • • • , 2g - 4, 2g - 2} (mod 2g) 

^4 = {2m |m = 1, • • • , (g - l)/2} = {2, 4, • • • , (g - 3), (g - 1)}, 
we have 

S = SiUS2US3U 84 = {m\m = !,■■■, 2Q}\{0,Q}. (mod 2g) (B7) 

Since P and 2g are mutually prime, S is invariant under a multiplication by P (mod 2g). 
This leads to 
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(Q-i)/2 



2m \ 



m=l 



n 

mg{l,-,2Q}\{0,Q} 



z — e 2Q 



[B8) 



In a limit z ^ 1, we have 



1^(1)1' = Q, 



(B9) 



which is independent of P. 

In order to obtain the other amplitude, we use a resursion relation which is obtained 
from Eq. {Em: 



1 _ q-2{i-i)-i I _ ^2(«-i)+i 



1 _ g-2(/-l)+Q-2 I _ ^2(1-1) 



1 _ g-2{/-l)-(Q+l) I _ ^2(Z-1)+{Q+1) 



1 _ g-2(Z-l)-2 I _ q2(l-l)+2 

sin7rg(2/ - 1) 



sin 7rT^2l 



From this we have 



[2/ - 1].!! 
[2/],!! 



(BIO) 



(Bll) 



where [2n]^\\ = [2] J4], ■ ■ ■ [2n], and [2n - 1],!! = [1],[3], ■ ■ ■ [2n - 1],. 

Now fix P = 1 and take a large Q limit. The continuum coordinate xi and the square 
amplitudes are defined by 



xi = il--)Ax, (/ = l,---,(Q + l)/2) 
2 



(B12) 
(B13) 
(B14) 



In the large Q limit, x is in an interval (0, 1). Taking Eq. ( |B10D upto first order in Ax and 
we get 



lognfx) 

ax 



d 

dx 



Sm TTX. 



(B15) 
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Thus 



n{x) 



C 



Sm TTX 



(B16) 



where C is a some constant. The coefficient C is determined by considering a case I 
0{1) « Q where 



fc-i 

n{xk) = Q n 

3=1 
k-l 



sin27rg(j 



sin27rQj 

^ • 1\ 2 
^-2] 



(BIT) 



APPENDIX C: DERIVATION OF THE FINITE SIZE CORRECTION FOR THE 

LARGE Q 

We derive the finite size correction near the edges of the wave function for P = 1 and 
E = 0. Write the ampfitude as {■^{xi)\^ = where X21+1 = {I + l/4)2/g. The 

finite size correction C{1) near the edge is given by 

X 

1 ,[2Z-1],!!,2 . 2tt.^ 1, 

- 2^^^^2>ll2k + l 4k' 
where / finite but is large but Q — > 00. 
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FIGURES 

FIG. 1. The energy spectrum of the Bethe Ansatz equations in the complex plane for (a) 
g = 23 and P = 1, 3, 5, • • • , Q - 3, Q - 1 (+) and (b) Q = 401 andP = 1, 3, 5, • • • , Q - 3, Q - 1 (• ) 

FIG. 2. The roots of the Bethe ansatz equation for the case P = 1 and Q = 89 in the complex 
plane. 

FIG. 3. The roots and pseudo roots of the Bethe ansatz equation for the rational fluxes which 
converges to 1/r, = Pk/Qk = 3/5,13/21,55/89,233/377, 987/1597,4181/6765 (a) in the whole 
complex plane and (b) an enlarged figure. In each case, the roots and the pseudo roots are always 
on the unit circle. The radii for 0^ are scalled to show the branching rule. 

FIG. 4. The roots of the Bethe ansatz equation for the rational fluxes which converges to 
1/(7, = Pk/Qk = 3/7,7/17,17/41,41/99,99/239,239/577, and 577/1393 (a) in the whole complex 
plane and (b) an enlarged figure. The radii are scaled. 

FIG. 5. Squared amplitudes for the wavefunctions for P = 1 and Q = 5,21,89,377 . The 
wavefunctions are normalized by the peak heights. 

FIG. 6. Amplitude of the wavefunctions: for the ratios of successive Fibonacci numbers: 
(a) whole region = 3/5, 13/21, 55/89 ) and (b) enlarged plot in the region near 0.89. 
{(f, = Pk/Qk = 3/5, 13/21, 55/89, 233/377, 987/1597, 4181/6765) The wavefunction with larger 
value of Q is shaded darker. The wavefunctions are normalized by the peak height. 

FIG. 7. /(a) for the wavefunction at = for the golden mean case. 

FIG. 8. The roots of the Bethe Ansatz equations for 7^ states in the complex plane for 
P = 1 and Q = 61. (a) The highest energy state, (b) the second higest energy state, (c) the third 
highest energy state, and (d) the 30-th energy state which is just one above the center {E = 0) 
state. 
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FIG. 9. All sets of the roots of the Bethe Ansatz equations in the complex plane for P = 1 and 
Q = 41. All the roots are on the unit circle. We scale the radii to show them together. The radius 
of the higher state scaled smaller. 

FIG. 10. Approximated distribution function of the roots of the Bethe Ansatz equations for 
P = 1 and Q = 4:1. (a) The highest energy state, (b) the 3-rd energy state, (c) the 10-th energy 
state, (d) the 15-th energy state, (e) the 18-th energy state, and (f) the 21-th energy state which 
is the center {E = 0) state. 

FIG. 11. The roots of the Bethe Ansatz equations for E states for P = 55 and Q = 89. (a) 
The highest energy state, (b) the 22-nd energy state, (c) the 30-th energy state, and (d) the 34-th 
energy state. 
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